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We calculate the rate of creation of double occupancies in a 3D Fermionic Mott insulator near half-filling by 
modulation of optical lattice potential. At high temperatures, incoherent holes lead to a broad response peaked 
at the Hubbard repulsion U. At low temperatures, antiferromagnetic order leads to a coherent peak for the 
hole along with broad features representing spin wave shake-off processes. This is manifested in the doublon 
creation rate as a sharp absorption edge and oscillations as a function of modulating frequency. Thus, modulation 
spectroscopy can be used as a probe of antiferromagnetic order and nature of quasiparticle excitations in the 
system. 



Advances in experiments with cold atoms on optical lattices 
have made them promising candidates for simulators of lattice 
models which play an important part in our understanding of 
strongly interacting quantum systems. Recently Mott insulat- 
ing statesJU have been obtained in the large U limit of the 
repulsive Fermionic Hubbard modelj2][3l with these systems. 

Motivated by recent experiments[2|, we give a theoretical 
formulation of the response of a Mott insulator near half- 
filling to modulation of optical lattice potential by relating the 
rate of production of double occupancies to the convolution of 
spectral function of holes and double occupancies (doublons) 
in the Mott insulator. We will show that this technique can be 
used to detect the presence of antiferromagnetic (AF) order 
and probe the nature of quasiparticle excitations (coherent vs. 
incoherent) in the system. We also discuss the connection of 
this response to optical conductivity in corresponding charged 
systems. 

We focus on two temperature regimes: (i) the high temper- 
ature limit ( J <C T r~j t h <c [/), where T is the temperature, 
th the tunneling matrix, and J = it^/U the super-exchange 
scale, which controls the quantum dynamics of the back- 
ground spins; and (ii) The low temperature limit (T <C J), 
with an AF ordered spin background. 

In the paramagnetic phase (current experiments!!]), we get 
a response peaked at to = U with a width equal to twice the 
bandwidth of the holes, reflecting the completely incoherent 
hole and doublon in this limit. We also derive a sum rule for 
the energy integrated rate of doublon production in this limit. 

At low temperatures, the AF ordering leads to coherent 
propagation of quasiparticles and manifests itself in a sharp 
absorption edge in the production rate. Additional structures 
at higher energies appear as a result of shake-off processes 
of spin waves. Thus, lattice modulation spectroscopy can be 
used to observe the nature of quasiparticle excitations and de- 
tect the presence of antiferromagnetic ordering in the Mott 
insulating state. 

Modulation of Optical Lattice: Repulsive fermions in op- 
tical lattices are well described by a one band Hubbard model 
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where the tunneling matrix th and the onsite repulsion U de- 
pend on the depth of the optical lattice V through |4| 
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where Er is the recoil energy of the photon, A is its wave- 
length and a s is the s-wave scattering length of the atoms. 

The modulation of the optical potential V(t) = Vq + 
SVsin(uit) effects both th and U. The modulation of U 
can be neglected in the large U /th limit. The modulation of 

th(t) =th + 5th sin(wt) is related to V(t) by 



6t h = t h 6V 
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The single particle spectrum of the Mott insulator is formed 
of two bands: (a) the lower Hubbard band, which does not 
contain double occupancies and is exactly filled at half-filling 
and (b) the upper Hubbard band, which contains a single dou- 
ble occupancy and is completely empty at half-filling. These 
are separated by the Mott gap ~ U. So, the modulation of 
the optical barrier will produce double occupancies once the 
frequency of modulation exceeds the Mott gap. 

Schwinger Bosons and Slave Fermions: In the Schwinger 
Boson Slave Fermion representation, we represent the singly 
occupied sites (spins) by two Schwinger bosons cft a , the dou- 
bly occupied sites by a doublon aft and the empty sites by 
a holon h\ The doublon and holon are Fermions. The 
original Fermion creation operator can be written as cJ CT = 
a\ a hi + cTai-c,d\ with the local constraint equation aft ia a; La + 
d\di + h\hi — 1. We define the following operators: FL = 

E CT 4r a i* and A lj = £ CT aa \v a )-<, which represent the 
hopping of the bosons and the creation of singlet configura- 
tions. We further make the following transformation: On B 
sublattice aft — > — eft . Then the unperturbed Hamiltonian is 

Ho = t h E(4 d j + h\ h i) F H + (4 h ] A v + h - c -) 
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FIG. 1 : (a) The density of states of a single hole in a half-filled back- 
ground in the atomic limit, (b) The rate of production of double oc- 
cupancies as a function of frequencies and (c) The energy integrated 
response in high temperature limit. 



while the perturbation due to the lattice modulation is 

H x it) = St h sm[tu t] 4 h ) A *3 + hdjAl (5) 

(H) 

where we have neglected terms that do not create or destroy 
doublons. We work with a system at half-filling, at tempera- 
tures T -C J7, where we can neglect the presence of doublons 
and holes in the unperturbed system. The number of doublons 
created in time t is given by 

N d {t) = J2^ En/T (n\UHt)J24diU(t)\n) (6) 

n i 

where \n) denotes the unperturbed states (in our case 
spin configurations) with energy E n and the time evolu- 
tion operator has a perturbation expansion U(t) = 1 — 

i /„* dt'H[(t') - f* dt' /*' dt"H[(t')H[(t"), where the stan- 
dard interaction representation of an operator is given by 
O 1 ^) = e lHot O(t)e~ iHot . We also assume that the dou- 
blons are created by the action of the perturbation Hamilto- 
nian only, i.e. the time evolution of the system by the unper- 
turbed Hamiltonian conserves the number of doublons. This 
neglects the decay of doublons into a pair of bosons during 
the time evolution. This approximation is justified as long as 
T <C U, since, due to energy conservation requirements, the 
decay of a doublon in the system is a very slow process 0. 
Under these assumptions, the first order response of the sys- 
tem vanishes and upto second order in perturbation theory, the 
rate of creation of doublons is given by 

P d {u) = \{5t h ) 2 J dwi J duj 2 Yu 

{ij){lm) 
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where uj is the frequency of the perturbation, A dt ^ is the 
spectral function for the doublon (hole) and A s is the Fourier 
transform of [Ai m (t)A\.(t')). It is to be noted that we have 
already used a mean field decoupling of the doublon, hole and 
Schwinger boson operators to arrive at the above equation. 
This approximation is justified at large U /th due to the sepa- 
ration of energy scales governing the hole (doublon) dynamics 
(th) and the spin dynamics (J). 

We emphasize that the response we are calculating is not 
equivalent to optical conductivity in the condensed matter sys- 
tems, (i) The current vertex in optical conductivity is re- 
placed by the kinetic energy vertex, (ii) The optical conductiv- 
ity involves convolution of hole spectral function with itself, 
whereas the calculated response involves convolution of hole 
and doublon spectral functions. Since the doublon spectral 
function is shifted by U, as we move away from half-filling, 
there is no response at low frequencies, whereas there would 
be optical response at low frequencies in a compressible state. 

High Temperature: We now focus on the regime U ^S> 
T ~ th ^> J, which is the regime of interest for the cur- 
rent experiments. In this limit, the quantum dynamics of the 
spins are irrelevant and one can replace the A operators by the 
probability of finding a TJ. or J.T configuration in an ensemble 
where all spin configurations occur with equal weight. Thus, 

P„( W ) = ~P s (St h ) 2 j du'Ai^A^-J) 

(ij) (Im) 

(8) 

where P s is the probability of finding relevant configurations 
at and (l,m) given by P s = 1/2 if = (l,m), 

P s = 1/8 if and (7, to) have no overlap and P s = 1/4 
if (i, j) and (I, to) share one site in common. Due to particle- 
hole symmetry of the problem at half-filing, we have A d {oj + 
U) = A h (uj), so that it is enough to compute the spectral 
function for the holes only. 

We now try to evaluate the spectral function of a single hole 
in a half-filled background where the spins are completely dis- 
ordered. The hole is completely incoherent, i.e. it moves dif- 
fusively in the system. The spectral function of the hole in this 
limit has been worked out by Brinkman and Rice 1 6 1 and Kane 
et al using the so called retraceable path approximation. 

The Green's function has contributions from processes 
where the hole hops from one point to another. However, as 
the hole hops, it scrambles the spin configuration and a string 
of ferromagnetic bonds is required along the path for the pro- 
cess to contribute. The probability of finding such a string is 
given by (1/2) L , where L is the length of the path. However, 
the trajectories where the hole retraces its path do not scram- 
ble the background spins and have a weight of 1 as opposed to 
(1/2) L . They provide the dominant contribution to the den- 
sity of states at low energies. 

The easiest way to derive the spectral function is to write 
the Green's function as a function of the frequency co in the 
following way : G -1 ^) = u>[l — E(u;)] and derive a self- 
consistent equation for the self energy. The first contribution 
to the self-energy comes from the hops to nearest neighbours 
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where z is the co-ordination 



and gives T,^\u>) 
number of the lattice. To include the longer hops, the denom- 
inator of E W can be modified with a higher order self-energy 
E( 2 )(w) = zt 2 h /u 2 (l - (z- l)t 2 /uj 2 ), where the factor z - 1 
comes from excluding the initial site while considering the 
initial hop. This method is similar to the one used by Ander- 
son in his original paper on localization physics 1 8 1 and gives 
a self energy E(w) = (z/z- l)[l/2- ^u> 2 - 4(z - l)t 2 /to}. 
This leads to the spectral function 
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The spectral function is plotted as a function of frequency in 
Fig[TJa). The spectrum is incoherent and has a band-width of 
2-s/z— lift. The spectral weight decreases monotonically as 
one goes towards the band edge. The rate of production of 
doublons, calculated using this spectral function is plotted in 
Fig |TJb). There is a peak around ui — U with weight upto 
twice the bandwidth (for the holes) around it. 

We note here some recent work in the paramagnetic phase 
using different techniques l9l [lOl . 

Sum Rule: Sum-rules have played an important role in var- 
ious strongly correlated systems, since they often involve less 
approximations and serve as useful check on both theory and 
experiments. In this case we consider the energy integrated 
rate of production of doublons J du>Pd(to). Using the identity 
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temperature limit 
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Sij we obtain the sum rule in the high 



duP d (u) = -{8t h ) 2 



The sum-rule is proportional to (8th) 2 , which is proportional 
to t\ for a constant fractional change in the amplitude of the 
lattice potential. Assuming U /t is tuned by tuning the lattice 
potential and a constant fractional change in the amplitude of 
the potential (SVo/Vo is held fixed), one finds that in the Mott 
regime, the energy integrated weight monotonically decreases 
with increase in Vq/Er, as shown in Fig |T|(c). 

Low Temperature (Antiferromagnetic phase) : We now 
consider the response of the system at T — in an AF ordered 
phase. As we will see, in this regime, the spectral function of 
holes has a sharp peak and a series of broad features. The 
sharp peak is reflected in the doublon creation rate as a sharp 
absorption edge and the broad features result in oscillations in 
the rate as a function of modulating frequency. 

This phase is characterized by Bose condensation of f and 
J, Schwinger Bosons on opposite (A and B) sublattices. In a 
1/5 expansion, the fluctuations are governed by the Holstein 
Primakoff Hamiltonian 



(11) 
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FIG. 2: (a)The spectral function of the hole at (it/2, n/2, n/2) (solid 
line) and at (0, 0, 0) (dashed line) for a system with U/t = 20. The 
spectral functions are broadened by an artificial broadening of 0.06t 
(b) The rate of doublon production Pd(ui) as a function of perturbing 
frequency. 



with 7k = (1/3) [cos k x + cos k y + cos k z ], and the quasipar- 
ticle operators are given by = — «k o^_ k with 



(12) 



(10) The coherence factors and i>k are given by = 



(l/v^K 1 + 1) 1/2 and v k = -sgn( 7k )(l/V2)(o 
l) 1 / 2 . The hole hopping Hamiltonian can be written as 

H h = t h ^2 h]h,a,„a' rT 

(ij) 



(13) 



Replacing the f and J. spins on A and B sublattice by the con- 
densate amplitude y/po = 1, this term can be written as 



H h = Zth^/PQ ^ fe k fek-q(Mq7k-q«q + WqTk^-q) + h.C. 
kq 

(14) 

The motion of a hole is thus accompanied by creation of a 
spin wave. We calculate the self-energy of the hole in a self- 
consistent Born Approximation ifTTI [121 which is equivalent 
to calculating the non-crossing Feynman diagram for the self 
energy. At T = 0, the self-energy is given by 

E(k, u) = |r(k,q)| 2 G(k - q, w - co q ) (15) 



where the vertex function T(k, q) = -rff lv / po(u q 7k-q 
v q 7k) and the self-consistency is ensured through 



where w k = z,J(l — 7 £) x l 2 is the dispersion of the spin wave 



G- X (k,u;) =w-53(k,w). 



(16) 
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The spectral weight obtained from the self-consistent solution 
for J = Q.2th is plotted for two different k values, ([0,0,0] 
and [it/2, it/2, tt/2]) in Fig^a). At the lowest energy of prop- 
agation of the hole, which occurs at (tt/2, tt/2, tt/2), spin 
waves cannot be created (at T = 0) and there is a coherent 
peak. The coherent weight is largest at this point and grad- 
ually decreases as one moves to the center of the Brillouin 
zone. The location of the coherent peak disperses as ~ J"fk> 
corresponding to second order hopping processes which do 
not scramble the AF alignment. 

Beyond the coherent peak, there are additional broad 
features at higher energies, whose peak to peak distance 
scales with J. These are generated by spin wave shake-off 
processes |[T3l . The peaks correspond to 2, 4, 6, ... spin waves 
and are dominated by spin waves near the Brillouin zone 
boundary where the fiat spin-wave spectrum results in a di- 
verging density of states. This is similar to peaks in the 2- 
magnon Raman response in antiferromagnetic insulators| 14|. 

In terms of the calculated spectral function, one can calcu- 
late the rate of doublon production as 

P d {w-U)= ^(5t h ) 2 K J2ll [ cL; 1 A(k,Lj 1 )A(k,ij-u 1 ) 

k 

(17) 

where k = 1 — (1 /2z) ^ k j^/cdk is a vertex correction which 
takes care of the singlet spectral function. There is no con- 
volution with spin spectral functions as the spin dynamics is 
governed by the transverse (phase) modes (f on B and J, on A 
sublattices) and the longitudinal (amplitude) modes (| on A 
and J, on B) of the condensate are neglected 0. 

The rate of doublon production is plotted is Fig |2|b). It 
shows an abrupt edge at the lower end of the spectrum cor- 
responding to the coherent spectral weight of the holes. The 
other oscillations in the response reflect the convolution of the 
coherent part with the broad incoherent peaks due to shake- 
off processes and that of the peaks themselves. We thus see 
that the presence of the AF order leaves its signature in the 
frequency dependence of the response. 

We sketch what happens as we move away from half-filling 
(still remaining within the AF phase). The order parame- 
ter po decreases, weakening the scattering of the hole by the 
spin waves. Thus the coherent part should grow leading to a 
sharper edge. This is in contrast to disordering the AF phase 
by raising temperature, where the scattering from occupied 
spin wave modes reduce the coherent part. 

Although our calculation is done for T = 0, we expect 
these qualitative features to be valid as long as the temperature 
is much below the Neel transition temperature. 

Comparison with Experiments: In the experiments of 
Ref. 12 J, the modulation is kept on for a fixed number of cy- 
cles. The time of drive is proportional to u^ 1 and the quantity 
measured is proportional to P^{iiS)/u). So the frequency in- 
tegrated response should decrease even faster with U/th as 
compared to Fig. [TJ;. However, the experimental data shows a 
monotonic increase with U/th\l5\. 

We believe this could be due to several reasons: (a) The 



response might be dominated by terms beyond second order 
perturbation theory. This can be checked by putting on the 
drive for different amount of time and looking at the linear- 
ity (or lack thereof) of the number of doublons produced with 
time, (b) The unperturbed system is not in thermal equilib- 
rium due to slow relaxation of the doublons created during 
tuning U/t Q or (c) Relaxation of doublons while driving 
the system leads to a steady state behaviour. We hope the dis- 
crepancies between the theory and experiments can be settled 
with further experiments on this system. 

Conclusion: We have related the rate of production of dou- 
ble occupancies by modulation of optical lattice in a 3D Mott 
insulator near half-filling to the convolution of hole and dou- 
blon spectral functions. This technique can be used to study 
the nature of quasiparticle excitations and detect presence of 
AF order in the system. In the paramagnetic phase there is a 
broad response peaked around u> = U. In the ordered phase, 
the coherent hole is reflected as a sharp absorption edge, while 
shake off processes lead to oscillations. 
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